Theory of Magnetic Field Induced Spin Density Wave in High Temperature 

Superconductors 
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The induction of spin density wave (SDW) and charge density wave (CDW) orderings in the mixed 
state of high T c superconductors (HTS) is investigated by using the self-consistent Bogoliubov-de 
Gennes equations based upon an effective model Hamiltonian with competing SDW and d-wave 
superconductivity interactions. For optimized doping sample, the modulation of the induced SDW 
and its associated CDW is determined by the vortex lattice and their patterns obey the four- 
fold symmetry. By deceasing doping level, both SDW and CDW show quasi-one dimensional like 
behavior, and the CDW has a period just half that of the SDW along one direction. From the 
calculation of the local density of states (LDOS), we found that the majority of the quasi-particles 
inside the vortex core are localized. All these results are consistent with several recent experiments 
on HTS. 



PACS numbers: 74.20.-z, 74.25Jb, 74.60. Ec 

Intensive efforts have been focused recently on the na- 
ture of vortex core excitations and the possible induc- 
tion of SDW and other phases in the mixed state of high 
T c superconductors (HTS). Experiments from the scan- 
ning tunneling microscope(STM) neutron scatter- 
ing |^-|| and nuclear magneticresonance(NMR) j?J pro- 
vided vital informations on these topics. For example, 
according to the neutron scattering experiment by Lake 
et al. @, a remarkable antiferromagnetism or SDW ap- 
pears in the optimally doped I^-aSr^CuC^ when a 
strong magnetic held is applied. Very recently, STM mea- 
surement by Hoffmanei al. || studied the LDOS in the 
mixed states of optimally doped Bi2Sr2CaCu20s+5, and 
they found that associated with SDW, an inhomogeneous 
CDW exist both inside and outside the vortex cores. The 
coexistence of d- wave superconductivity (DSC) and SDW 
orders was theoretically studied by several groups in the 
absence of a mag netic field §-[□]]. In the presence of 
a magnetic field, this problem was studied by the SO (5) 
theory and also by a Ginzburg-Landau approach ]l3| . 
The LDOS at the vortex core in a pure DSC were first 
calculated by authors in Ref.(14). With induced anti- 
ferromagnetic (AF) order, the LDOS was investigated 
without the contribution of the quasiparticles 
and by a mean field study [jl5| . 

Although the competition between the SDW and DSC 
in a magnetic field was previously examined, the nature 
of the induced SDW and its spatial variation have not 
been addressed in such detail as to compare with the 
experiments. In this paper, we shall adopt the method 
described in previous papers Jl4|,|l5| to examine the pos- 
sible induction of extended SDW and CDW orders in the 
mixed state of HTS, and their nature in under-, optimally 
and overdoped samples. In order to simplify the numer- 
ical calculation, we shall assume a square vortex lattice 
for the mixed state and a strong magnetic field B such 
that A 3> b ^> £ with A as the London penetration depth, 
£ the coherence length and b the vortex lattice constant. 
Under this condition, the applied magnetic field B can be 



regarded as a constant throughout the sample. Our cal- 
culation is based upon an effective model Hamiltonian 
with competing SDW and DSC orderings. In the ab- 
sence of a magnetic field B and for the optimally (hole) 
doped sample x = 0.15, the parameters are chosen in 
such a way that the SDW ordering is completely sup- 
pressed and only the DSC ordering prevails. When the 
system is in the mixed state driven by a magnetic field B, 
an inhomogeneous SDW is induced. We found that the 
structure of the induced SDW is determined and pinned 
by the underlying vortex lattice. For optimized doping 
sample, the modulation of the induced SDW and its as- 
sociated CDW is determined by the vortex lattice and 
their patterns obey the four-fold symmetry. With the 
same set parameters and B=0, SDW, DSC and CDW 
stripes appear in the underdoped sample (x — 0.12). In 
the presence of a strong B, all order parameters pattern 
show remarkable anisotropic behavior along the x and y 
directions, and the CDW has a period just half that of 
the SDW along one direction. Increasing the SDW in- 
teraction strength can also lead to quasi-one dimensional 
pattern. The LDOS near the vortex core has also been 
calculated. It shows an asymmetric double peaks around 
E = in agreement with the results of Ref.(15) and ex- 
periments QJ^] . In addition our results also show a small 
gap at E = indicating the presence of the SDW order 
in the bulk sample. From the spatial profile of the LDOS, 
we conclude that almost all the quasiparticles inside the 
vortex core are localized. 

Let us begin with a phenomenological model in which 
interactions describing both DSC and SDW order param- 
eters in a two-dimensional lattice are considered. The 
effective Hamiltonian can be written as: 

+ J2(A^4 r cl l + h.c). (1) 
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The summation here is over the nearest neighbors sites. 

c\ a is the electron creation operator and /i is the chemical 

potential. In the presence of magnetic field B, the hop- 

ii i /*' A ( r )' dr 

ping integral can be expressed as ty = toe r -> 

for the nearest neighboring sites For simplicity, 

we have set the next-nearest neighbor hopping equal to 
zero. $o = h/2e is the superconducting flux quanta. 
Here we choose Landau gauge A = (—By, 0,0) with y 
as the y-component of the position vector r. The two 
possible orders in cuprates are SDW and DSC which 
have the following definitions respectively: Af' DW — 
J7(4 T c iT - CjjjCil) and Ay = V DS c{c\\c^ - c u c jT )/2. 
In the above expressions, U and Vjjsc are respectively 
the interaction strengths for SDW and DSC orders. The 
mean- field Hamiltonian (1) can be diagonalized by solv- 
ing the resulting Bogoliubov-de Gennes equations self- 
consistently 
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where the single particle Hamiltonian Tifi = — ty + 
(Uriia — //)<5ij , and 

ni ,=J2K\ 2 f(E n ), (3) 
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with f(E) as the Fermi distribution function and the 
electron density n\ = m-\ + rtij . The DSC order parame- 
ter at each site i is Af = (AP +Bx i + AP_ ex i - Af i+Gy - 

Afi_ ey )/4 where APj = A-^exp^ j^ )/2 A{v) -dr] 
and e^ ^ denotes the unit vector along (x,y) direction. 
The main procedure of self-consistent calculation is given 
below: For a given initial set of parameters n- la and Ay, 
the Hamiltonian is numerically diagonalized and the elec- 
tron wave functions obtained are used to calculate the 
new parameters for the next iteration step. The calcu- 
lation is repeated until the relative difference of order 
parameter between two consecutive iteration step is less 
than 1CP 4 . The solutions corresponding to various dop- 
ing concentrations can be obtained by varying the chem- 
ical potential. 

In the following calculation, the length and energy are 
measured in units of the lattice constant a and the hop- 
ping integral to respectively. We need to point out that 
the induction of internal magnetic field by the supercur- 
rent around the vortex core is so small comparing with 
the external magnetic field that we can safely adopt the 



unfiorm magnetic field distribution approximation. We 
follow the standard procedures |l4|,[l5| to introduce mag- 
netic unit cells, where each unit cell accommodates two 
superconducting flux quanta. By introducing the quasi- 
momentum of the magnetic Bloch state, we obtain the 
wave function under the periodic boundary condition 
whose region covers many unit cells. The related pa- 
rameters are chosen as the following: The DSC coupling 
strength is Vdsc = 1-2, the linear dimension of the unit 
cell of the vortex lattice is chosen as N x x N y = 40 x 20 
sites and the number of the unit cells M x x M y = 20 x 40. 
This choice corresponds the magnetic field B ~ 37T. 

At optimal doping and B = 0, the important qual- 
itative characteristic of the quasiparticle states is iden- 
tical to that of a pure d-wave superconductor. At the 
first step, the magnitudes of parameters are selected to 
fulfill such requirement that the AF order is completely 
suppressed and only DSC prevails in the absence of a 
magnetic field. Here we choose the hole doping x = 0.15 
and U — 2.4. Our calculation is performed at very low 
temperature where the vortex structure is almost inde- 
pendent of temperature. In Fig. 1(a) we plot the typical 
configuration of vortex structure. The DSC order param- 
eter vanishes at the vortex core center and recovers its 
bulk value at a couple of coherence lengths away from the 
center. By comparing it with the vortex structure of a 
pure DSC, the size of the vortex core here is noticed to be 
enlarged. The centers of the two vortex cores are at sites 
(10,10) and (30,10). Fig. 1(b) displays the spatial dis- 
tribution of the staggered magnetization of the induced 
SDW order as defined by M? = (-l)*Af DW /U. It is 
obvious that the AF order exists both inside and out- 
side the vortex cores, and behaves like a two-dimensional 
SDW with the same wave length in the x and y directions. 
The induced SDW order reaches its maximum value at 
the vortex core center and the magnitude of its spatial 
variation still holds the fourfold symmetry as that of the 
pure DSC case. The order of DSC and SDW coexist 
throughout the whole sample. The appearance of the 
SDW order around the vortex cores strongly affects the 
spatial profile of the local electron density distribution, 
which can be represented by a weak CDW as shown in 
Fig. 1(c). The remarkable enhancement of electron den- 
sity (or depletion of the hole density) is presented at the 
vortex core center. It is easy to observe that the four-fold 
symmetry holds for both SDW and CDW with the same 
wavelength 20a. 

To have a deeper understanding of the above results, 
the case of an underdoped sample [x — 0.12) is exam- 
ined. After the calculation is performed at B = 0, we 
found that DSC, SDW and CDW orderings have the 
stripe structuress, which is consistent with stripe phase 
results reported by Martin et al. |ll[] . We have compared 
the free energy between stripe phase solution with a uni- 
form AF solution and found that the free energy of the 
stripe phase is always lower. In the presence of a strong 
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magnetic field B, the profile of DSC order parameter is 
shown in Fig. 2(a). The radius of vortex is further en- 
larged than in Fig. 1(a). The qualitative features of Fig. 
2 are quite different from those in Fig.l and the spatial 
variation of the SDW becomes quasi-one dimensional. Its 
periods of oscillations are fixed by the vortex lattice sim- 
ilar to the case of x = 0.15. We also noticed that the AF 
order is much enhanced in the mixed state as compared 
with its values in the stripe phase at B = 0, in agreement 
with Katano's experiment 0. The results are presented 
in Fig. 2(b) and Fig. 2(c), where the anisotropy in the 
magnitudes of SDW and CDW along the x and y direc- 
tions shows quasi-one-dimensional behavior. The periods 
of SDW and CDW are respectively 20a and 10a. These 
results are in qualitative agreement with the observations 
of Lake et al. ^ where a magnetic field induced striped 
AF order was observed in underdoped sample. We also 
have calculated the case for overdoped sample (x = 0.20) 
with the same set of parameters. For B = 0, the SDW 
order is completely suppressed and the DSC order is ho- 
mogeneous in real space. When B is strong, the SDW 
order does not show up even inside the vortex core. 

We notice that the periods of SDW and CDW obtained 
from the experiments are respectively 8a and 4a for 
optimal doping sample. With the present band param- 
eters, we are not able to obtain these numbers. How- 
ever, the experimental values could be obtained by tuning 
U, doping or including a next-nearest neighbor hopping 
term. For larger U case, the configurations of SDW and 
CDW exhibit stripe-like behavior along the x— (or y— ) 
direction. The periodicity of CDW is always half that 
of SDW. For example when U = 3 and x = 0.20, the 
periodicity of SDW and CDW are respectively 10a and 
5a. But under this condition, the AF stripe phase would 
appear in the optimally doped sample. 

Next we present the calculation of the LDOS near the 
center of vortex core for hole doping x = 0.15. The LDOS 
is given by: 

Pi (E) = -Y,M\ 2 f'(E n — E) + \v?\ 2 f'(E n + E)] , (6) 

n 

It also measures the differential tunnel conductance, 
which could be measured by STM experiments. In this 
case, the thermally broadening effect has not been con- 
sidered here and the temperature is fixed at T = 0.01. 
We plot the LDOS at the core center vs the quasi-particle 
energy measured from the Fermi level. For comparison, 
we have also displayed the LDOS at the site between 
two next-nearest neighboring vortex cores. The results 
are shown in Fig. 3. At the core center, the coher- 
ent peaks due to the gap edges of the superconductor 
at B = are suppressed, and two asymmetric peaks of 
the vortex states appear slightly above and below E = 0. 
The spectrum agrees qualitatively with the experiment 
for YBCO and that of Ref.(15). But a closer inspec- 
tion reveals that the distance between these two peaks 



becomes slightly larger than that in Ref.(15). The peak 
at about E = 0.8 seems to be the characteristic of the 
bulk SDW order and the band structure effect. It is easy 
to notice that the magnitude of the LDOS at E = ap- 
proaches zero which indicates a bulk SDW gap. From the 
spatial distribution of LDOS at E = (see Fig. 4(a)), we 
find that LDOS is enhanced along the x — y and x = —y 
or the diagonal directions from the center of the vortex 
core. This 'star'-like behavior indicates the small number 
of quasiparticles very close to E = are extended [jl6| . 
We also present the result (see Fig. 4(b)) for the peak 
on the left of E = in Fig. 3. Its profile displays an 
obvious localized shape even though it decays somewhat 
slower along the x and y directions than the diagonal 
directions. The peak on the right at E ~ 0.16 has a sim- 
ilar spatial distribution which is not shown here. From 
Fig. 4(b), we conclude that the majority of quasiparti- 
cles inside the vortex core with energies near or at the 
two asymmetric peaks are in fact localized, in agreement 
with experiments . 

In conclusion, we have studied the induction of SDW 
and CDW orders in optimally and underdoped HTS by 
a strong magnetic field. Consider only the nearest neigh- 
bor hopping term, the spatial variations of DSC, SDW 
and CDW orders have been numerically presented in Fig. 
1 and Fig. 2 and stripe-like structure exhibits for under- 
doped sample. The LDOS near the vortex center have 
also been calculated. We show that almost all the quasi- 
particles inside the core are localized, that is very differ- 
ent from case for a pure DSC vortex These results 
are consistent with recent STM, neutron scattering and 
NMR experiments on HTS. Finally we would like to em- 
phasize that although our self-consistent BdG calculation 
based upon a mean-field approach tends to overestimate 
the stability of the SDW phase, the qualitative features 
of our results should still be valid, particularly in view of 
their favorable comparisons with experiments. 
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FIG. 1. The amplitude distribution of the DSC order pa- 
rameter AP (a), the staggered magnetization M? (b), and 
the electron density ri\ (c) in one magnetic unit cell. The 
size of the cell is 40 x 20, corresponding to a magnetic field 
H = $o/(20 x 20). The strength of the on-site repulsion 
U — 2.4 and the averaged electron density n = 0.85. 



FIG. 2. The amplitude of the DSC order parameter Af 
(a), the staggered magnetization Mj (b), and the electron 
density n\ (c) in one magnetic unit cell. The averaged electron 
density n — 0.88. The other parameter values are the same 



as Fig. 1. 



FIG. 3. Quasiparticle LDOS profiles at the center of the 
vortex core are shown by solid lines. The dashed lines show 
the LDOS at the site midway between two next-nearest neigh- 
bor vortices. The averaged electron density n = 0.85. 



FIG. 4. The spatial distribution of LDOS in a vortex unit 
cell at the energy £ = 0a) and E = -0.126 (b). 
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Fig. 1(a) Chen & Ting 




Fig. 1(b) Chen & Ting 




Fig. 1(c) Chen & Ting 




Fig.2(a) Chen & Ting 




Fig.2(b) Chen & Ting 




Fig.2(c) Chen & Ting 




Fig. 3 Chen & Ting 
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Fig. 4(a) Chen & Ting 




Fig.4(b) Chen & Ting 




